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Abstract 

In this paper, the real scalar field equation in Schwarzschild-de Sitter 
spacetime is solved numerically with high precision. A method called 
'polynomial' approximation is introduced to derive the relation between 
the tortoise coordinate x and the radius r. This method is different from 
the 'tangent' approximation [1] and leads to more accurate result. The 
Nariai black hole is then discussed in details. We find that the wave 
function is harmonic only near the horizons as I. Brevik and B. Simonsen 
[1] found. Howerver the wave function is not harmonic in the region of 
the potential peak, with amplitude increasing instead. Furthermore, we 
also find that, when cosmological constant decreases, the potential peak 
increases, and the maximum wave amplitude increases. 

1 Introduction 



I. Brevik and B. Simonsen [1] solved the real scalar field equation in Schwarzchild- 
de Sitter spacetime numerically. They approximated the tortoise coordinate 

x = x(r) by the 'tangent' function r = r(x) and replaced the potential function 

v = v(r) by v = v(r) — v(x). The tortoise coordinate x was introduced to 
simply the real scalar field equation. Then they solved the equation by Mat- 
lab software, and found that the solution in SDS system is close to that of a 
harmonic wave [1] . 
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In the present paper, we continue the previous work and investigate the 
solution of the real scalar field equation in Schwarzchild-de Sitter spacetime. 
In section 2, we will introduce a new method, referred to as the 'polynomial' 
approximation, to get the relation between the tortoise coordinate x and the 
radius r. This approximation is different from the 'tangent' approximation 
introduced by I. Brevik and B. Simonsen [1], and approximates the curve of the 
tortoise coordinate x versus the radius r more accurately. Then we solve the 
real scalar field equation by Matlab software. The solution illustrated in figures 
demonstrates that the solution of the real scalar field equation in Schwarzschild- 
de Sitter spacetime is not harmonic globally in contrast to the result from the 
work in Ref.(l). 

We adopt the signature ( — h ++), put h, c, and G (Newton's gravitational 
constant) equal to unity, and follow the conventions of Misner et al.[2] 

2 Global Solution of the Real Scalar Field Equa- 
tion in Schwarzschild- de Sitter spacetime 

2.1 Horizons 

With our conventions Einstein's equations read 

R ^ - \ R 9»v + A 9^ = 8ttT^. (1) 
For a spherically symmetric system the line element takes the form 

ds 2 = -f(r)dt 2 + --^-dr 2 + r 2 dB 2 + r 2 sin 2 6d(p 2 . (2) 
f{r) 



In our case 



f(r) = 1 — , (3) 

r 3 

where M is the mass of the black hole. I. Brevik and B. Simonsen[l] have 
discussed that when AM 2 < |, equation /(r) = has three real solutions. They 
are event horizon r e , cosmological horizon r c and a negative solution r a , with 
r e > 0, r c > 0, r Q = -(r c + r e ). The choice AM 2 = 0.11 leads to r e = 2.8391M 
and r c = 3.1878M. 

2.2 The Wave Function 



The real scalar field equaion in this case can be written as 



= - 



,tt 



f(r) 



+ 



1 

r 2 sin 9 



(sin6»$0) ; 



sm( 



0. (4) 
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Let 



$ = -y u i(r)e- i " t Y lm (e,<p), 
r 

where Yi m is the usual spherical harmonic. We have 



v(r) 



(5) 



(6) 



where v(r) is the potential 
v ( r ) = f(r) 



f'(r) 1(1 + 1) 



/ 2M Ar 2 \ (2M 2A 2 \ , N 



^ is taken as unity here. We now introduce the tortoise coordinate x by the 
equation [3] 



1 f dr 



x = 



2Mj f( r y 



(8) 



This quantity is conveniently expressed in terms of the surface gravities k^, 
defined by [4] 
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We get 
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A, 



A. 



(9) 

(10) 

(11) 
(12) 



The tortoise coordinate can now be written as 



2M 



J-lnf^-! 



J-Infl-il 



(13) 



From Eq. (8), we get 



dx/dr = 



2/W 



(14) 
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So Eq. (6) is rewritten as 



•* Lul (x)=4MWV^(x), (15) 
where v is now a function of x. 

v = v(x). (16) 
2.3 Boundary Conditions and Results 

To proceed further, we need the inversion r = r(x) of the tortoise coordinate 
x = x(r) to get v = v(r) = v(r(x)) = v(x). But it is difficult to invert Eq. 
(13). In order to get the function r = r(x), we turn to the numerical way. For 
Eq. (13) is a one to one mapping function, we use the command 'polyfit' in 
Matlab software to get a new function y = y{x) = Y]"—n ciiX 1 , n G N, i G N, to 
approximate tortoise coordinate r = r(x). In different interval of the radius r, 
the function y has different form. In the same interval, n is determined by the 
demand of our approximation accuracy. The bigger n is, the higher the accuracy 
is. Now we put M = 1, A = 0.11 and consider Nariai black hole [1]. Thus the 
horizons arc r e = 2.8391 and r c = 3.1878. In the interval r G [2.83908,3.18775], 
we have 

20 

y = ^2a lX \ (17) 

i=0 

{cii} are displayed in Table 1. 

The curve of y versus the tortoise coordinate x and the curve of the radius r 
versus the tortoise coordinate x are illustrated in Figure 1, from which we know 
that the two curves overlap with each other in the interval r G [2.83908, 3.18775] . 
In figures, point A denotes the position at which the potential peak is. Now 
we see that the function y may be better as an approximation to Eq. (13) 

than the tangent function [1] r — 15tan[6(r — d) + 5], with b = 2.7 /(r c — r e ), 

d = (r c + r e )/2. The curve of r versus the tortoise coordinate x and the curve 
of the radius r versus the radius x were illustrated in Rcf. (1). 

Now inserting y = y{x) into form (7) instead of r = r(x), we obtain the 
potential v as a function of x 

v = v(x). (18) 

The curve of v versus x from equation (13) and the curve of v versus x from the 
approximation (17) is illustrated in Figure 2. The curve of v versus y and the 
curve of v versus the radius r are illustrated in Figure 3. Figure 2 and Figure 3 
show that our approximation is good. 

We set v = near the horizons. Thus Eq. (15) becomes 

* w jO»0 = 44Ux), (19) 



dx 2 



+ AM z v 




4 



which has the following solution 



®wi(x) = cos(2a;). 



(20) 



Therefore, we have [1] 



100) = = 100) = cos(200). 



(21) 



The above numerical values are used in the boundary conditions when Eq. 
(15) is solved with the Matlab software. The solution is illustrated in Figure 4 
and Figure 5. These two figures show that the wave function near the horizons 
is harmonic, but it isn't in the region where the potential gets its peak value. 
There is a largest amplitude at x = 3.1431 near the point x = 2.8582. While 
the potential peak v = 7.5 x 10~ 4 appears at x = 2.8582. When r trends to the 
point x = 3.1431, the amplitude of the wave increases step by step. I. Brevik 
and B. Simonscn [1] illustrated the partial solution of Eq. (15) and mentioned 
that the solution is close to that of a harmonic wave. Now we find that the wave 
function is not harmonic globally. 

The case when A = 0.001 will not be discussed here again. Its result is 
similar to the one when A = 0.11. But one thing should be mentioned. It 
is that, When A decreases, the potential peak increases. So does the largest 
amplitude of the wave. 



In this paper, we introduce a new method called 'polynomial' approximation to 
approximate the tortoise coordinate x = x(r). And we find that the solution of 
the real scalar field equation in Schwarzchild-de Sitter spacetimc is harmonic 
near the horizons as I. Brevik and B. Simonsen [1] found, and that it isn't near 
the position where the potential gets its peak value, with higher wave peaks and 
lower troughs appearing instead. 

3.2 About the Reflection and Transmission Coefficient 

In this paper, we did not calculate the reflection and transmission coefficients. 
The reason is that there are many square potentials [1] to describe a common 
potential v — v{x). Each square potential gives out a reflection coefficient and a 
transmission coefficient by the method adopted in I. Brevik and B. Simonsen's 
paper [1] . Different square potential gives out different reflection coefficient and 
transmission coefficient. 
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Summary 



3.1 



What We Have Done 
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Tablel. the coefficients in function y = y(x) 



a = 2.9817 


ai = 6.5107 x 10" 3 


a 2 = 4.0912 x 10"° 


a 3 = -2.9913 x 10~ b 


a 4 = -3.4895 x 10" 8 


a 5 = 1.6009 x 10" y 


a 6 = 2.3413 x lO" 11 


a 7 = -8.0083 x 10~ i3 


a 8 = -1.1964 x 10~ 14 


a 9 = 3.3845 x 10~ 16 


aio = 4.3110 x 10~ 18 


an = -1.0899 x 10" 19 


a w = -1.0120 x 10" 2i 


ai 3 = 2.4364 x 10" 23 


aw = 1.4031 x 10" 25 


ai 5 = -3.4329 x 10" 2 ' 


ai 6 = -9.5242 x 10 _3U 


ai 7 = 2.6439 x 10~ 3i 


ai 8 = 1.5652 x 10~ 34 


ai 9 = -8.0402 x 10" 3e 


a 20 = 4.3 2 6 2 x 10" 3y 
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Figure 1: The curve of the tortoise coordinate x versus r (full line), together 
with the curve of the tortoise coordinate x versus polynomial approximation y 
(dotted). A = 0.11. 
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Figure 2: The curve of the potential v versus the tortoise coordinate x (dot- 
ted), together with the curve of the potential v versus x from the polynomial 
approximation y = y(x) (full line). A = 0.11 
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Figure 3: The curve of the potential v versus the radius r (dotted), together 
with the curve of the potential v versus y (full line). A = 0.11. 
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Figure 4: The curve of the wave function ty^i versus x in SDS. A = 0.11. 
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Figure 5: The curve of the wave function $> u l versus r in SDS. A = 0.11. 
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